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Squark and gaugino hadroproduction and decays 
non-minimal flavour violating supersymmetry 
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Abstract. We implement non-minimal flavour violation (NMFV) in the MSSM at low scale. In 
this framework, we propose benchmark points for the mSUGRA scenario including non-minimal 
flavour violation by evaluating a number of experimental low-energy, electroweak precision and 
cosmological constraints. We finally discuss phenomenological aspects of our NMFV scenario and 
present a numerical analysis of squark and gaugino production cross sections at the LHC. 

PACS. 12.60.Jv Supersymmetric models 



1 Flavour violation in the MSSM 



While the only source of flavour violation in the 
Standard Model (SM) arises through the rotation of 
the quark interaction eigenstates into the basis of phys- 
ical mass eigenstates, in Supersymmetry (SUSY) ad- 
ditional flavour violation is introduced at the weak 
scale through renormalization group running pQ. In 
minimal flavour violating (MFV) scenarios, the non- 
diagonal entries of the squared squark mass matri- 
ces stem from the trilinear Yukawa couplings of the 
fermion and Higgs supermultiplets and the resulting 
different renormalizations of quarks and squarks. In 
constrained MFV (cMFV) scenarios, these flavour vi- 
olating entries are neglected at the SUSY breaking and 
weak scale. 

However, when SUSY is embedded in larger struc- 
tures such as GUT theories, new sources of flavour 
violation arise [2], that have to be included in the 
squared squark mass matrices at the weak scale. As 
the corresponding flavour violating off-diagonal entries 

AfJ cannot be simply deduced from the CKM-matrix 
alone, in NMFV is parametrized by considering them 
as free parameters. Denoting Mj- , M R and rrikXk 
their usual diagonal and helicity mixing entries, k = 
1,2,3 referring to the (s)quark family, the squared 
squark mass matrices are then given by 
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The scaling of the flavour violating entries A qq 
with the SUSY-breaking scale Msusy implies a hier- 
archy Af' L > A qq ' LLR > A qq R [2J. The numerically 

largest A qq L are also related by SU{2) gauge invari- 
ance through the CKM-matrix, not allowing a large 
difference between them. They are usually normalized 
to the diagonal entries [3] , 



A qq 



Kj M lq M Jq 



(2) 



so that NMFV is gouverned by 24 arbitrary dimen- 
sionless parameters X qq . 

The diagonalization of the mass matrices Mj and 
Mj requires the introduction of two additional 6x6 
matrices R u and R d , which relate the helicity eigen- 
states to the physical mass eigenstates through 

(ui, . . . ,u 6 ) T = R u (u L ,c L ,t L ,u Rl c Rl t R ) T ', (3) 
0i,..., d 6 ) T = R d (d L , s L ,b L ,d R , S R , b R ) T . (4) 

By convention, the squark mass eigenstates are la- 
belled according to < ... < rriq e for q = u,d. All 
relevant couplings are then generalized by expressing 
them in terms of R u and R d . 



2 Constraints on NMFV in mSUGRA 

In order to take implicitly into account a number of 
experimental bounds coming from the neutral kaon 
sector, B- and D-meson oscillations, rare decays, and 
electric dipole moments [SJH], we restrict ourselves to 
the case of flavour violation between only the sec- 
ond and third generations of left-chiral squarks, imple- 
mented through one real NMFV-parameter A = \ 8 £ L = 
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\ c l L , while all other XfJ are zero. Allowed regions for 
this parameter are then obtained by imposing explic- 
itly low-energy, electroweak precision and cosmological 
constraints. We start with the theoretically robust in- 
clusive branching ratio [7] 

BR(6-> s 7 ) = (3.55 ±0.26) x 1CT 4 , (5) 

which affects directly the allowed squark mixing bet- 
ween the second and third generation. Since NMFV 
contributes already at one-loop level, as do the SM 
contributions, this constraint is supposed to depend 
strongly on the NMFV-parameter A. A second observ- 
able, sensitive to squark mass splitting and thus to 
NMFV, is the electroweak p-parameter, for which new 
physics contributions are constrained to [1] 

Ap = -uT = (1.02 ±0.86) x 1(T 3 . (6) 

by the latest combined precision measurements. A third 
variable sensitive to SUSY loop corrections is the ano- 
malous magnetic moment of the muon, for which we 
require the new physics contribution a^ USY to close 
the gap 

Aa^ = (22 ± 10) x 10~ 10 (7) 

between recent experimental data and the SM predic- 
tion [3] . Finally, to have a viable cold dark matter can- 
didate, we require the lightest SUSY particle (LSP) to 
be neutral in electric charge and colour, and the re- 
sulting neutralino relic density ^cdm^ 2 to agree with 
the observational limit [12] 

0.094 < fi CDM /i 2 < 0.136. (8) 

The above limits are imposed at 2a confidence level 
on the mSUGRA model with five free parameters mo, 
mi/2, Ao, tan/3, and sgn(/i) at the GUT scale. The 
renormalization group equations are solved numeri- 
cally to two-loop order using the computer programme 
SPheno 2.2.3 [13]. At the weak scale, we generalize the 
squark mass matrices by introducing the NMFV pa- 
rameter A as described above, diagonalize these ma- 
trices and compute the low-energy and electroweak 
precision constraints using the computer programme 
FeynHiggs 2.5.1 [14]. The neutralino relic density is 
computed using a NMFV-adapted version of the com- 
puter programme DarkSUSY 4.1 |llj . taking into ac- 
count the six-dimensional hclicity and flavour squark 
mixing. For the Standard Model input parameters, i.e. 
the masses and widths of the electroweak gauge bosons 
and quarks, the elements of the CKM-matrix, the SM 
CP-violating phase, and Fermi's coupling constant, we 
refer the reader to Ref. [4]. 

The one-loop SUSY contributions to a^ USY are pro- 
portional to sgn(p) [15], so that the disagreement be- 
tween experiment and SM prediction is increased for 
/i < in all SUSY models. As furthermore, this region 
is virtually excluded by the b — > sj constraint for all 
A G [0;0.1] except for very high SUSY masses, we re- 
strict ourselves to the case of p > 0. The dependence 
on Ao of the variables is extremely weak, so that we 
use only Aq = throughout our analysis. 



Table 1. Benchmark points for mSUGRA allowing for 
flavour violation among the second and third genera- 
tions. We also indicate the allowed range for the NMFV- 
parameter A. The values of mo, in 1( / 2 and Ao are given in 
GeV. 
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Fig. [T] shows a typical scan of the mSUGRA pa- 
rameter space in mo and mi/2 for tan/3 = 10, Ao=0, 
and /i > 0, and different values of A G [0;0.1]. As ex- 
pected, the b — > S7 excluded region depends strongly 
on flavour mixing, while the regions favoured by a M 
and the dark matter relic density show quite low sen- 
sitivity on the A-parameter. Ap constrains the param- 
eter space only for very high scalar (mo > 2000 GeV) 
and gaugino (m 1 / 2 > 1500 GeV) masses, so that the 
corresponding regions are not shown. 

Based on an extensive analysis of the mSUGRA 
parameter space at different values of tan/3 (10, 30, 
and 50), we define benchmark scenarios, which are al- 
lowed/favoured by the above constraints, permit non- 
minimal flavour violation among left-chiral squarks of 
the second and third generation up to A < 0.1, and 
are at the same time "collider-friendly" , i.e. have rel- 
atively low values of too and m\/2- Our choices are 
presented in Tab. [TJ together with the correspond- 
ing allowed ranges for the NMFV-parameter A. Note 
that all four benchmark points are also valid in MFV 
(A < 0.005 . . . 0.01) and cMFV (A = 0) scenarios. 



3 Numerical discussion 

In this talk, we focus on the phenomenology of bench- 
mark scenario B. For this point all sparticle masses 
are below 1 TeV (see Fig. [3]), which makes it the most 
"collider-friendly" among the four proposed scenarios, 
having light sleptons (mj ~ 200... 300 GeV), relatively 
light gauginos (m-^ ~ 150. ..550 GeV) and squarks (m^ 
~ 650... 850 GeV), as well as a heavy gluino (m^ ~ 900 
GeV). For a similar discussion of the other points pre- 
sented in Tab.[TJ the reader is referred to Ref. [To] . 

The explicitly imposed constraints are shown in 
Fig. [2] except for the SUSY contribution to the anoma- 
lous magnetic moment of the muon, a^ USY ~ 14 x 
10 -10 , that is independent of A and agrees with the 
limits in Eq. ([7]) . The first graph shows the electroweak 
precision observable Ap. On our logarithmic scale, only 
the upper bound of the 2a range is visible. Despite the 
strong dependence of Zip on squark flavours, helicities 
and masses, the experimental error is sufficiently im- 
portant to allow for relatively large values of A < 0.52. 
The most stringent low-energy constraint comes from 
the b — > sj branching ratio, which depends strongly 
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Fig. 1. The (m , m 1/2 )-planes for tan/3 = 10, A = 0, [i > 0, and A = 0,0.03,0.05,0.1. We show a M (grey) and WMAP 
(black) favoured as well as 6 — > S7 (blue) and charged LSP (beige) excluded regions of the mSUGRA parameter space 
in constrained minimal (A = 0) and non-minimal (A > 0) flavour violation. 



Ap l> -> » T Oqa, n 2 LSP 




Fig. 2. Dependence of the precision variables BR(6 — > sj), Ap, and the cold dark matter relic density QcDMh 2 as well 
as the lightest neutralino, up- and down-type squark masses on the NMFV-parameter A in our benchmark scenario B. 
The experimentally allowed ranges (within 2<r) are indicated by horizontal dashed lines. 



on the NMFV-parameter A due to the same loop-level 
of the SUSY and SM contributions. The small exper- 
imental error allows for two narrow intervals in A. As 
the second one is disfavoured by B — > X s [i + [i~ [17], we 
restrict ourselves to the region A < 0.1. The neutralino 
relic density ^cdm^ 2 agrees with the limits in Eq. JSJ 
up to high values of A < 0.9, which is due to the fact 
that in the relevant neutralino (co)annihilation pro- 
cesses squarks only appear as propagators or at one- 
loop level. The decreasing of .!?cdm^ 2 for higher A is 
due to the larger mass splitting of the squarks, im- 
plying that the lightest up- and down-type squarks 
become lighter and coannihilation processes become 
more important. 

Fig. [3] shows the mass eigenvalues of the up- and 
down-type squarks as a function of the NMFV param- 
eter A <G [0;0.1]. We observe large mass splitting for 
the lightest and heaviest squarks with increasing A due 
the growing relative importance of the off-diagonal ma- 
trix elements. However, the masses of the intermediate 
squarks are practically insensitive to A. Another inter- 
esting phenomenon observed is the so-called "avoided 
crossing" between the squark mass eigenvalues. For 
the down-type squarks for example, this is observed 
at A w 0.02 between the states d\ and d 3 (see Fig. [3]). 
At the point where two levels should cross, the corre- 
sponding squark eigenstates mix and change character 
(see Fig. 0} . The same phenomenon occurs also for the 
up-type squarks at A ~ 0.02 between H4 and M5, and 
a second time at A ~ 0.1 between U3 and U4. These 
"avoided crossings" , a common phenomenon in quan- 
tum mechanics, are due to the fact that the Hermitian 
squark mass matrix depends continuously on one sin- 
gle real parameter A. 
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Fig. 3. Dependence of the up- and down-type squark 
masses on the NMFV-parameter A in our benchmark sce- 
nario B. 



In Fig. 2] we show the hclicity and flavour decom- 
position of some of the up- and down-type left-handed 
squark mass eigenstates for the experimentally favou- 
red range in the vicinity of (c)MFV, i.e. A 6 [0; 0.1]. 

We finally provide numerical predictions for pro- 
duction cross sections of squarks and gauginos at the 
LHC, i.e. for pp-collisions at y/s = 14 TeV centre-of- 
momentum energy. Fig. [5] shows the numerical results 
for some interesting channels in our benchmark sce- 
nario B. Again, for a complete discussion of all chan- 
nels in all four benchmark scenarios, the reader is re- 
ferred to Ref. 16J. The magnitudes of the cross sec- 
tions vary from barely visible level (10~ 2 fb) for weak 
production of heavy final states over the semi-strong 
production of squarks and gauginos and quark-gluon 
final states to large cross sections (10 2 to 10 3 fb) for 
the strong production of diagonal squark- (anti)squark 
pairs or weak production of very light gaugino pairs. 
Unfortunately, since the strong gauge interaction is in- 
sensitive to (s)quark flavours and gaugino pair produc- 
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Fig. 4. Dependence of the chirality (L, R) and flavour (u, c, t; d, s, and b) content of the discussed up- and down-type 
squark mass eigenstates on the NMFV-parameter A £ [0; 0.1] for benchmark scenario B. 




Fig. 5. Cross section examples for charged and neutral squark-antisquark pair, associated squark-neutralino, and 
gaugino-pair production at the LHC as function of the NMFV-parameter A € [0;0.1] in our benchmark scenario B. 



tion cross sections are summed over exchanged squark 
flavours, the processes with the largest cross sections 
are practically insensitive to the flavour violation pa- 
rameter A. 

Some of the subleading, non-diagonal cross sec- 
tions, however, show sharp transitions. For example, 
at A = 0.02, the cross sections for didg and d^d^ switch 
places, corresponding to an avoided crossing in the 
mass eigenvalues and a switch in the flavour contents 
of the di and d^ squarks (see Figs. [3] and 0J. Since 
the concerned sstrange and sbottom mass differences 
are rather small, this is mainly due to the different 
strange and bottom quark densities in the proton. At 
A = 0.035, the cross sections for d^d^ and d\d\ increase 
sharply, corresponding to the flavour content change in 
the g?3 squark from sr to da (see Fig. 01, that can now 
be produced from down valence quarks. The d\d^ pro- 
duction cross section increases with the strange squark 
content of d\ . 

Smooth transitions and semi-strong cross sections 
of about 1 fb are observed for the associated produc- 
tion of third-generation squarks with charginos and 
neutralinos. The cross section for d§ production de- 
creases with its strange squark content, while the bot- 
tom squark content increases at the same time (see Fig. 
0}. The opposite happens for g?3, where the produc- 
tion cross section increases. For the up-type squarks, 
the cross section of uq production increases/decreases 
with its charm/top squark content. 

4 Conclusion 

We have presented an implementation of NMFV in 
the MSSM at the weak scale. The relevant flavour 
violating off-diagonal elements of the squared squark 



mass matrix are parametrized by one real NMFV pa- 
rameter. Based on an extensive analysis of low-energy, 
electroweak precision and cosmological constraints, we 
proposed benchmark points for the mSUGRA scenario 
including NMFV. We discussed the phenomenology 
of these benchmark scenarios, which are characterized 
by the phenomenon of "avoided crossings" . We finally 
presented numerical predictions for squark and gaug- 
ino production cross sections at the LHC, among which 
several processes are sensitive to NMFV. 
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